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The equations of transient stimulated Raman scattering on the finite interval are solved by 
the spectral transform method on the semi-line. As the problem has a free end, the pump and 
Stokes output at finite distance can be constructed as the solution of a hnear Cauchy-Green integral 
equation. 



I. INTRODUCTION 



^ . 

Since the first observation in 1983 of the creation of a Raman spike in the pump depletion zone by Driihl, Wenzel 
' and Carlsten , there has been a large amount of studies of stimulated Raman scattering (SRS) of long laser pulses in 
<U , gas. The reason is that the Raman spike has been shown to be the macroscopic manifestation of large fluctuations of 
' the phase of the initial Stokes wave. Raman spike generation were then predicted have been given a coherent-mode 
^sj description , and experiments have been preformed where SRS grows spontaneously on initial fluctuations Q . The 
04 Raman spike hence appears as a means to study the quantum properties of the Stokes wave initiation, which gives 
informations on the phase of the electromagnetic vacuum. This comes in addition to the process of Stokes growth 
which amplifies the quantum fluctuations of the medium. A general discussion of the quantum coherence properties 
of SRS is given in g. 

As the SRS equations possess a Lax pair , there has been many attemps to modelize the Raman spike as a soliton 
^) ' (see e.g. ||]). But it has been proved recently that actually the Raman spike occurs in the spectral transform 
, theory as a manifestation of the continuous spectrum (hence it is not a soliton) when for a short period of time the 
reflection coefficient becomes close to zero 0. These results were obtained by solving the initial-boundary value 
' problem for the SRS equations on the infinite line, which, from a physical point of view, is inconsistent with a finite 
dephasing time of the medium oscillators. However the results obtained are strikingly close to the experimental data 
^d| . Here we construct a solution of SRS on the finite interval by using the spectral transform on the semi-line, first 
proposed in the context of nonlinear polarization dynamics . 

The interaction of light with a material medium, in the case when a laser pump pulse (frequency lul, envelope 
Al) interacts with the optical phonons (eigenfrequency uy, envelope Q) to give rise to a down-shifted laser pulse 
O (Stokes emission, frequency ujs, envelope As) according to the selection rules (Aw is the the detuning from the Raman 
resonance, Ak is the phase missmatch) 



UJs = LOl ^ ^ ^"-^j kg ^ ki^ — ky — Afc (I"l) 



. can be modelized by the following slowly varying envelope approximation jisf 

(4 + -4)^^ ^ ^^V^JiZJI QAs exp[z(AfcZ - Ac.T)], 
oZ col 4 rjc 

+ -4)^5 = ^^V^JiZJZ Q* Al exp[-z(AfcZ - AujT)], (1-2) 
oZ col 4, r]c 

4fQ + ^Q^ '-^.[^ AlA*s exp[-z(AfcZ - Ac^T)]. 
ol J 2 4 mujv V 

The electromagnetic wave field E{Z,T) and the material excitation X{Z,T) are obtained as 

E{Z, T) = ]-Al exp[i(fciZ - wlT)] + ]:\[^As exp[z(/cs^ - ^sT)] + c.c. (1-3) 
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X{Z, T) = exp[i(fcyZ - LovT)] + c.c. (1-4) 

Hereabove is the differential polarizability at equilibrium, c/rj is the light velocity in the medium, N is the density 
of oscillators of mass to, T2 is the relaxation time of the medium. 
For a medium initially in the ground state 

Q(Z,0) = 0, (1-5) 
and for an arbitrary set of input pump and Stokes envelopes profiles (for any value of the missmatch Afc) 

Al{Q,T) = Il(T-'^Z), As{Q,T)^Is{T-'^Z), (1-6) 
c c 

we give here the output values of both light waves profiles Al{L,T) and As{L,T) in terms of the solution of a 
Cauchy-Green linear integral equation. 

The solution, constructed by the inverse spectral transform theory (1ST), is actually exact for infinite relaxation 
times, and we have proposed an approximate solution for finite T2 which matches the experiments with high accuracy 
[ pT| . The method applies also for a non- vanishing initial state of the medium but the solution in this case is in general 
not explicit (this is relevant in physical situations where the quantum fluctuations of the initial population density of 
the two- level medium is taken into account). 

The fact that 1ST can be applied to SRS on the finite interval has been first proposed by Kaup in 1983 jsj. However 



the evolution of the spectral data given there does not correspond to the boundary problem {1-6) and in particular, 
as this evolution is homogeneous, it does not allow for the growth of the Stokes seed on a medium initially at rest. In 
a different context, a nonhomogeneous evolution of the spectral data has been obtained in where the self- induced 
transparency equations are solved for an arbitrary initial (t = —00) population density of the two-level medium. 1ST 
has been later used to solve an initial-boundary value problem on the half-line for the nonlinear Schrodinger equation 
(NLS) by Fokas and Its iQ, but in this case the required boundary data in x = for the potential itself renders 
nonlinear the evolution of the spectral transform. 

The property of SRS of being solvable on the finite interval results simply from the nature of the equations for 



which the initial-boundary value problem (I-£)(I-6) is well posed and does not require new constraints when it is 
given on the finite interval (this is not so for NLS for which the vanishing boundary values at infinity become some 
prescribed boundary value in a; = 0). Consequently the method applies for every other case of solvable evolutions 
with nonanalytic dispersion relations when precisely passing to the finite interval does not imply adding information 
or constraint. However we will discover that there are requirements of analyticity of the boundary data in order for 
the problem to be integrable. 



Before going to the method of solution, it is convenient to rescale the system (1-2) into a dimensionless system, 
which goes with defining the new variables 

x^jZ, t=j{T-^Z). (1-7) 
L L c 

Then the dimensionless rescaled material excitation is defined as (the differential polarizability has the dimension 
of a surface) 

qix, t)^- ^^^JZ^ LQ{Z, T), (1-8) 
4 77c 

while Al ans As are rescaled by using the boundary conditions as 

AL{Ak,Z,T) As{Ak,Z,T) 
aL[k,x,t)^ , as(k,x,t)^ , /„ = max|/L|, (1-9) 

which scales to 1 the incident laser pump amplitude. The dependence on the phase missmatch is represented here by 
the dimensionless wave number 

fc = ^L(Afc- ^Aw), (I-IO) 
which we shall refer to as the essential missmatch parameter. Finally, for 

„ Na'pep ujs 2 ... 
16 ijmc^ ojv 
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(dimensionless) , the system ([-2) becomes for infinite dephasing time 



dxCtL = qas e 



-iA.U)t 2ikx 



dxcts = -q*aL e 



(1-12) 



„ „ „ „* ,jil^u]t —2ikx 

qt — —goaLas s s 



(1-13) 



Due to the dependence of the field envelopes and 05 on the essential missmatch k, it is necessary to consider 
the cooperative interaction of all fc-components with the medi um. Then we need actually to consider the medium 
excitation as resulting from all k components and replace (H3) with 



qt = -.90 / dkaLttg e e 



(1-14) 



where now the input [x — 0) values of the pump wave and the Stokes wave as are also function of k sharply 
distributed around fc = (the resonnance), that we denote by 



aL{k,Q,t) = JL{k,t), as{k,0,t) = Js(fc,t). 



(1-15) 



We will demonstrate that the system ( I-12| )( I-14 ), with the initial data q{x,0) in and boundary values ( [-15 ), 
where JL(k,t) (resp. Js(k,t)) has an analytic continuation in Im(fc) > (resp. Im(fc) < 0) vanishing as |fc| 00, is 
integrable on the semi-line x > Q. Moreover the solution furnishes the output values of the pump and Stokes waves 
va. Z = L (i.e. a; = 1), which gives the solution of the SRS equations on the finite interval. 



II. THE SPECTRAL PROBLEM ON THE SEMI-LINE 

Basic definitions. We briefly give the basic notions on the Zakharov-Shabat spectral problem on the semi-line for 
the 2x2 matrix u{k^x^t) in the potentials q{x,t) and r{x,t) 



I/, 2: > 0, t> 0. 



(II-l) 



The solution v can be verified to obey 



dx 



det{v} = 0. 



Two fundamental solutions, say ly , are defined by 



^21 



Jo dx'qv+^ 



{11-2) 



(II-3) 



^12 

'42 



'12 



(II-4) 



Jo dx'qv^^ 



dx'rv^^e 



2ik{x — x') 



(II-5) 



f^dx'qv22'i ^''"'■'^ 



(II- 
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Riemann- Hubert problem. In order to derive a well-posed Riemann-Hilbert problem for i^(fc), we require that iy{k) 
be bounded for large k and hence define it as 



. _ ( i^+{k,x,t), Im(fc) > 
> ~ \iy-{k,x,t), lm(fc)<0 



(II-7) 



Indeed it obeys 



fc — > oo 



Kfc)~i + o(^), 



(II- 



which is obtained from the integral equations by integration by part. The scattering coefficients (functions of k and 
parametrized by t) are defined for k real, as 



/■oo 

/ dx' q4^e'''^' 
Jo 



p 



poo 

/ dx' 



(II-9) 



T+ = 1 + / dx' rv^^, = 1 



dx' qv2i- 



(11-10) 



Following standard methods (see e.g. especially the appendix), one can prove from the above integral equations 
the the solution v obeys on the real axis, the following Riemann-Hilbert problem 



(II-U) 



'1 ■ 

The method consists simply in comparing e.g. the integral equations for {y^ — i/^)e^^'^''^ to the one for ■ 

Bound s tates . As the integrals run on the finite support [0,a:;], the solutions and of t he Volterra integral 
equations (II-3)( |l-^ ) are analytic (entire functions of fc) in the complex plane. The property (11-2) is used to compute 
the determinants of [v^ ,V2) and of {^i tI'2) both at x = and x = 00, which gives 



I'liik, 00, t) 



1 



T+(fc,t)' 



, h'22{k,CO,t) 



1 



r-(fc,t)- 



(11-12) 



Hence the quantity l/r+ (l/r ) is an entire functions of fc and can have a number iV+ {N ) of zeroes fc+ (fc„ ) which, 
for bounded r and q are simple an of finite number. 



Consequently, the solutions and 1^2 of the Fredholm integral equations (II-5) and ([1-4), which can be written 
also 



ir dx'qi^2l 
dx'rv^-^e^'^^^^^^"' 



(11-13) 



^12 
^22 



,+ _(Cdx'q^t2e-^^^^^-^') 
ir dx'rv+2 



(H-14) 



have respectively the TV and simple poles fc„ and fc+ of t and r+ . Then the integral equations for the residues 
allow to obtain 



Res = C„ V2 (fc„ ) exp[2ifc„ x], C„ = Res p 

fen 



(H-15) 



Res U2 — C^i^J*"(fc^) exp[— 2ifc+a;], — Res 

fc+ fc+ 



(H-16) 



As for the Riemann-Hilbert relations (II-ll), the method consists in comparing the integral equations for 
exp[— 2ifc~a;]Res and for C'~t'^(fc~). 
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Boundary behaviors. An important relation for the following is the so-called unitarity relation 

T^T^ = 1 — p~ , 



obtained by computing i'22{k,oo,t) by using (11-11) and by comparing the result with ( |II-12| ). 
Consequently the boundary behaviors of ly can be written 



X = 



1 p+ 
1 



1 
p- 1 



(11-17) 



(11-18) 



,+ - ( .Vr+ 



" ' l/r- 



The d problem. The analytic properties of v can be summarized in the following formula 



with the spectral transform 



dv 

— V i? exp[2ifc(T3x], 

ok 







where the distributions and S{k — kn) are defined as 



j:c+s{k-k^„ 





(11-19) 



(11-20) 



(11-21) 



dX A d\fiX)S^iXi)^ -2i I dXnJiXR ± zO), 



dX A dXf{X)S{X - kn)= .f{kn) 



with the notation A = Xr + iXi. 



III. INVERSE SPECTRAL PROBLEM 



(11-22) 



The Cau chy- Green integral equation. The inverse problem is solved by integrating the 9-equation ( II-20D with the 
boundary (II-8). We prove here the following theorem: the solution f{k,x,t) of the Cauchy-Green integral equation 



fik,x,t)^l + 



2i7r 



dX A dX 

X-k 



/(A, X, t)R{X, t) exp[2iA(T3a::] 



(III-l) 



coincides with the solution v{k^x^t) of ( |II-3| )-(II-6) if the reflection coefhcient possess a meromorphic continuation 
with simple poles k^ and residues C^, namely 



d 



±lMk)>0 ^ ^p^ik)^2i7:J2ct Sik-kt: 



(111-2) 



Such analytical property of the spectral data occur when the potentials are on compact support. Therefore we shall 
assume here that the initial datum q{x, 0) is indeed on compact support and we will have to demonstrate that the time 
evolution conserves this analytical property. Actually the physically interesting case is when ^(x, 0) vanish identically. 
This implies no bound states which will be assumed in the following (everything can be extended to the case with 
bound states as they appear as the poles A;+ of p"*" in Im(fc) > 0, and those (fc~) of p^ in Im(fc) < 0). 
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Proof of th e Theorem. The first step is to verify from (0-20) that the function / is solution of the differential 
problem ([W). This is easily done by following the method of |15|, in short: prove the relation 



d_ 

dk 



[(/. - ikfas)/-'] = 



and integrate it. Then the solution / of ( |III-l| ) solves the spectral problem (|Tlj) with 

(° ?)=■[-/'■'] 

where Z*-^^ is the coefficient of 1/fc in the Laurent expansion of f{k). 
The second step consists in proving 

f^{k,x,t) = iy^{k,x,t), f2{k,x,t) ^ iy2{k,x,t), 



(III-3) 



(III-4) 



(III-5) 



which is acheived just by comparing the values of these vectors in a; = 0. The functions /"'' and / solve the following 
coulped vectorial system (forget for a while the {x, i)-dependence) on the real axis Im(fc) — 0: 



/i+(fc) 



dX 



1 



2i7r 7 A - (fc + iO) 



2i\x 



(III-6) 



f2ik) 



dX 



1 



2iT: J A - (fc - iO) 



— p+(A)/+(A)e-2'^- 



(III-7) 



Now, by closing the con tour o f integration in Im(A) < for and Im(A) > for /2 , the Cauchy theorem for the 
analyticity requirement (111-2) leads to 



Va; < , Im(fc) = : /+(fc, x, t) = 



(fc, x, t) — 



(III- 



Consequently, the two m atrice s (/]^,/2~) and {'^^^,'^2) solve the same first order differential problem and have the 
same values in x = 0, so (III-5) is proved. 
The third step results in proving 



fi {k,x,t) = {k,x,t), f^{k,x,t) = iy^{k,x,t), 
which is performed by evaluating the functions 



(III-9) 



1 



— r,- 



dX 



+ f+ -2i\x 



2m j A - (fc + «0) 
dX 



•^21 2iTT J A - (fc - iO) 



p-f22e^'^' 



as X — *■ c». A usefuU formula is 



lim P 



dX 



s'(^-fe)-/(A) = ±mf{k), 



x^±(X3 J X — k 

where Pj denotes the Cauchy principal value integral, and the Sokhotski-Plenielj formula 

dX 



A - {k±iO) 



f{X) = ±z7Tf{k) + pJ ^/(A). 



(III-IO) 



(III-ll) 



(III-12) 



Since it has been already proved that — and /2 = V2 ^ the functions /j^ and f22 have from ([L4) and ( |ll-5| ) 
bounded behaviors as x ^ 00. Consequently 



1™ /12 = 0. 1™ /21 = 0- 



(III- 13) 
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The remaining functions f^^ and are obtained by using the differential form of (III-l), namaly 



fn fii 



P /l2i /22 /; 



22 



-2ikx + r + 
P J2l- 



With ([Il-g), this gives 



Vx < , Im(fc) = : f^^{k,x,t) = 1, f+{k,x,t) = 1. 



(III-14) 



(III-15) 



FinaUy, the two matrices {fi,f2) ^^'^ {'^it'^'z) solve the same first order differential problem and have the same 
values in X = for the diagonal elements and in x = oo for the off-diagonal, hence they are equal. The theorem is 
proved and an immediate consequence is that the coefficients are effectively given by (|lI-9[). 



From now on we will use the single matrix as denoting the common solution to ( [ff-f| ) and ([I-2)-(II-£) 
Reduction. It is convenient for physical applications to assume the following reduction 



for which 



y =iyf2, (l^l^2 ) * = " ^^21 : 



(III- 16) 



(III- 17) 



The overbar denotes the complex conjugation and the "star" 



nk) = f{k). 

As a consequence the boundary values of v become 



X = : iy~^ 



1 P 
1 



1 

-p* 1 



(III- 18) 



(III-19) 



(III-20) 



X = oo : V 



l/T 







l/r* 



(III-21) 



where we use the notation p = and r = for real k. 

Transmission coefficient. It is usefuU for the following to derive the relationship between r and p. Although this 
can be done in general, we stil l cons ider only the case when no bound states are present (this will be the case of 
interest). From the definitions ( 11-10 ), we have for large k 



T^{k)^l+0{y) 



(III-22) 



Defining then h{k) a s ln(T +) for Im(fc) > 0, and — ln(r ) for Im(fc) < 0, its discontinuity on the real fc-axis can be 
written by means of ( 11-17 ) 



In(l-pV)- 



(III-23) 



As from (III-22), h{k) vanishes for large k, the above Riemann-Hilbert problem has the following solution 



lm(fc)^0 : h{k) = ^ j ^ ln(l-p+p-), 



and hence r are given from p . In the reduction (fII-16), this relation can be written 



d\ 



271 J A - fc 



-PI 



ln(l 



(III-24) 



(III-25) 
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IV. LAX PAIR AND GENERAL EVOLUTION 



Lax pair. The compatibility 

Ut-V^ + [U, V]=0 (IV-1) 

between the two following spectral problems 

ly^^iy ikaz + Uv, U ^ -ikas + ( ^ q ) , (IV-2) 



leads to the evolutions 



dt 



R^[R , n] + M. 



(IV-5) 



This general result is not proved h ere (s ee [p^ or ||l6|l ), but we give in the next section the computation of t he tim e 
evolution of p^{k, t) (given here by ([V-5)) directly from the Lax pair, when the evolution of q{x, t) is given by (IV-10). 

Solvable evolution. The entries n{k,t) (diagonal matrix-valued function) and M{k,t) (off-diagonal matrix-valued 
distribution) are arbitrary, which allows to solve evolutions with nonanalytic dispersion relations {dfl/dk ^ 0) and 
arbitrary boundary values (the problem of arbitr ary b oundary values has been first solved for the SIT equations 
in 1^). As we are int erested in solving evolutions as (1-14) where the integral runs on the real fc-axis, we work within 
the reduction (111-16) and chose (remember the definition ( [n-ig| ): 



M{k,t) = 
and the dispersion relation 

n{k,t) =uj{k,t)a3, uj{k,t) 



m{k,t)5^ 
m*{k,t)6'- 



dX 



(p{\,t), k^K, 



analytic except on the real axis where it stands the discontinuity 



^4r^ = i - i^-{kR,t)) 5{ki) - ^{kR,t)5{ki). 

ak 2 ^ 



(IV-6) 



(lV-7) 



(IV-8) 



Then the quantity v{d^/d\)v ^ is ill-defined because v{k, x, t) itself is discontinuous on the real axis. To compute it 
we need to use the identity Jl^ 



ok ok 



(lV-9) 



Within the reduction and with the above choice of M and il, the evolution ( [V-4 ), becomes by means of ( lV-9 ) 

dk [2iputi^i2 + m{i^ti?e~^''''' + m*(i/f2)^e2*'=^] . (IV-10) 



To get the above result, it is necessary to use the Riem ann-H ilbert relations (11-11) to express i^ii and i>{2 in terms 
of ly^i and z/j^. This equation, coupled to the system ( IV-2| ), is the general solvable e volut i on w ith a nonanalytic 
dispersion law vanishing at large k. It is used now to solve the boundary value problem ( t-12 )( l-14 ). 
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V. SOLUTION OF SRS 



The above theory aUows to solve the nonhnear evolu tion problem ( I-12| )([-14) with the initial data of q{x, 0) (which 
actually will be taken to vanish) and boundary value (|l-15| ). This is done in the following way 
Basic relations. As the 3 functions 



as exp[2ikx — iAcot] 



^21 



■^22 



exp[2i/ca;], 



(V-1) 



solve the sam e first-order differential system, they are uniquely related by their values in a; = and comparing ( I-15| ) 
with ( |[II-20| ) we obtain readily 



as eicp[2ikx — iAujt] 



= Jl 



^21 



Js 



^22 



exp[2iA;x — i/S.Ljt]. 



(V-2) 



This formula actually gives the solution to the physical problem (compute the output values from the input values) 
as soon as the function p^{k, t) is calculated (t is expressed from p in (III-25)). 

Computation of p{k, t). This computation is performed by first determining the functions i^ik, t) and m{k^ t). This 
is done by equating the evolut ion (1-14) with ( IV-10| ) by means of t he r elation (V-2) between {a^^as) and v. Using 
the reduction relations (III-17), we fin ally get that the evolution (1-14), where (a^, as) is expressed in terms of v 
through ( |V-2| ) , reads exactly as ( tV-lC ) for 



-^^9o{\Jl\'-\Js\ 



The corresponding evolution (IV-5), which reads 

p+ = -2oj^ 



2im, 



then gives by means of (IV-T) 

pt{k,t) = TrgoJL{k,t)Js{k,t)e 



'^-'-^90P^{k,t) 



dX 



X-{k + iO) 



i\jL\'-\Jsn 



(V-3) 



(V-4) 



(V-5) 



>5o 



dX 



X - kr, 



{\Jl?-\Js?) 



(V-6) 



For everything to work, we have seen in sec. 3 that the analytical requirement (III-2) is essential. It is clear on the 



above time evolution that this is ensured for all t if we require that JL{k,t) be analytic in the upper half plane and 
Js{k, t) in the lower one (note from (III-19) that Jg{k, t) is a function of k analytic in the upper half plane). In short 



Im(fc) > 



d_ 

dk 



J^(fc,t) = 0, Im(fc) < 



d_ 

dk 



Js{k,t) = 0. 



(V-7) 



The other constraint is that, for all t, the function p'^{k,t) vanish at large k in Im(A:) > 0. This is guaranteed here 
because JL{k,t), Js{k,t) and uj~^{k,t) do vanish at large k in Im(fc) > 0. Note that it would not be true for an 
analytic dispersion relation (like the polynomial k^ for NLS). 

A medium initially at rest corresponds to q{x,0) — 0, and hence from ( |ll-g| ) and ( 11-16 ) to 



p+(fc,0) = 0, C+(0) = 0. 



(V-8) 



Consequently the evolution ( |V-q ) ensures that no bound stat es are create d and the equation ( V-5 ) can be solved 
explicitely. Then the solution of the related integral equation (III-6) (III-7) gives u^{k,x,t) which in turn allows to 
compute the light field envelopes aL{k, x, t) and as (A:, a;, t). In particular their values in x = 1 furnishes the output in 
the case of the finite interval of physical length L. For instance, with the solution p'^{k,t) of (V-5), the pump output 
is obtained by solving (with p^ = —{p'^)*)- 



^ti{k,x,t) = 1 



1 

2i7r 



dX 

A - (fc + iO) 



p (A, t)z^j^2('^i 2;, i)e 



2i\x 



(V-9) 
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2i7r J 



dX 



P^(A,i)i^n(A,x,t)e' 



A - (fc - iO) 
and then it is given by 

aL{k, l,t) = JL{k,t)i^^^{k, l,t) + Js{k,t)i^^2{k, l,t)exp[2i/c - lAwi]. 



(V-10) 



(V-11) 



We note finally that, in the case of the semi-infinite line {L oo), from the bomidary values ( III-21 ), we have the 
following pump output 



(fc, oo, t) = — ( Jl - /9+ J; 



(V-12) 



which has been used in |pj| to interpret the experiments of 



VI. EVOLUTION OF THE SPECTRAL TRANSFORM FROM THE LAX PAIR 



The auxiliary Lax operator given in (IV-3) can be simplified by using ([V-9) and the property that lyflv ^ vanishes 
as 0{l/k). Then (|V^) reduces to 



ZITT 



dX A dX 

X-k 



^ — iXa^x 



{M+[R, n])e 



For the particular choices ( IV-6 ) and ( IV-7 ), the above equation reads 

1 r dX 



Vt ~ X jy, X{k, X, t) 



X-k 



X{\x,t), 



(VI-1) 



(VI-2) 



X{k,x,t) = (m - 



^21''21 '^11''21 



''iVn 1 cxp[-2ifca;] 



+ (m* - 



^21^12 '^12'^22 



\2^\2 
V7r^ v7 



exp[2i/ca::]. 



(VI-3) 



Our purpose here is to recover now the time evolution ( V-4 ) of the reflection coefficient p{k,t) by means of the 
usual method which consists in evaluating ( VI-2| ) at one boundar y, say in x = 0. The bound ary v alues of both 
and X~ in x = are easily calculated from those of and ly^ in ( |II-18| ). Then the equation ( |VI-2[ ) gives for and 
u~ successively 



Pt = 



dX 



TT J X- (k + iO) 



{iuj p — to). 



(VI-4) 



TT J X~ {k + iO) 



i-iuj (p+)*+to). 



(VI-5) 



dX 



TT J X- (k-iO) 



{iu!^ p'^ — m), 



(VI-6) 



dX 



1 



IT J X- (k-iO) 



— (-icj (p+)*+to). 



(VI-7) 



First, these equations are compatible if = (w^)*, whi ch is guara nteed through ( IV-7| ) by cp £ iR (see ( V-3 )). 
Then the above four equations reduce to two and we select (VI-4) and (VI-6) whose solution goes in two steps. 
First step: the equation ( VI-6| ) is realized if the function itu'^p'^ — m is analytic in Im(fc) > and vanishes as fc — > oo. 
This is true by construction for ijj{k,t) and this is true also for m{k,t) because we have required that JL{k,t) be 
analytic in the upper half plane and Js{k,t) in the lower one (rmember that from (III-19) that Jg{k,t) is a function 
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of k a nalytic in the upper half plane). Finally this holds also for the reflection coefficient p^(k,t) by the requirement 

(HI- „ „ . ^ 

Second step: by sub strac tion of (VI-4) and (VI-6), and by use of the Sokhotski formula ( III-12| ) we arrive precisely at 
the time evolution ( V-4 ) . 

Finally, using the behaviors in x = oo instead of a; = 0, furnishes the time evolution of the transmission coefficient 
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